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Abstract 

The infrared behaviour of the n-point functions of a Yang-Mills theory with 
a charged scalar field in the fundamental representation of SU(N) is studied 
in the formalism of Dyson-Schwinger equations. Assuming a stable skeleton 
TiJ- expansion solutions in form of power laws for the Green functions are obtained. 

^ For a massless scalar field the uniform hmit is sufficient to describe the infrared 

^ scaling behaviour of vertices. Not taking into account a possible Higgs-phase 

it turns out that kinematic singularities play an important role for the scaling 
O solutions of massive scalars. On a qualitative level scalar Yang-Mills theory 

yields similar scaling solutions as recently obtained for QCD. 
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1. Introduction 



As we believe today there are four fundamental forces in nature, namely the 
gravitational, the electromagnetic, the weak and the strong force. One promis- 
ing approach to unification of the different interactions is to describe the elec- 
tromagnetic, weak and strong force by means of quantum field theory, which is 
usually called the standard model of elementary particle physics. Formulating 
gravity as a quantum field theory is not possible the usual way, but since the 
effects of gravity are much weaker on short distances compared to the ones of 
the other three forces the standard model provides a good description for high 
energy physics. 

Within the standard model the gauge theory which specifies the strong inter- 
actions is called quantum chromodynamics (QCD), involving the dynamics of 
quarks and gluons. For high energies the theory is asymptotically free, as shown 
in the nobel prize work by Politzer, Gross and Wilczek [1]. That means, that 
the coupling of the quarks and gluons vanishes for infinitely high momenta. 
Therefore perturbation theory, which assumes a free theory being disturbed by 
a small interaction, is justified for high energies. 

As the QCD /3-function shows, the strength of the coupling increases with de- 
creasing momenta, and therefore perturbation theory cannot be applied any- 
more in this region. Thus there is a need for non-perturbative methods to find 
the origin of typical phenomena, as e.g. dynamical chiral symmetry breaking 
or confinement, which means the absence of colour-charged particles in asymp- 
totic states. Until today both features are not fully understood, neither their 
origin nor their possible relation to each other. For a review on the confinement 
problem and several possible explanations see [2]. 

A promising approach to the whole energy regime is the formalism of Dyson- 
Schwinger equations (DSE), see [3|, HI |5] for the original works or [6l [3 [H [9] 
for reviews about the application of DSEs in hadron physics. Within this 
method the Landau gauge turned out to be the "simplest" gauge, because of 
the transversality of the gluon propagator and the non-enhancement of the 
ghost-gluon vertex [ini [H] in the infrared. The formalism of DSEs in Landau 
gauge has been very successful in the study of the infrared properties of both 
the Yang-Mills sector and QCD [T2l[l3l[Tll[l5l[l6l[I7l[l8l[T9l[20l|2llE^ 

[2ai23[2nil2Zl[2Hl[2SlEniElE21E31ISl[Ml[^ 

Also other tools prove useful in the analysis of Landau gauge QCD, as e.g. 
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functional renormalisation group equations (FRGs) pT j l39 | HT j H2| Il3 | li H H5 | 
llHlllTlllHlliglEOllSTlEalESllMll^ Gribov-Zwanziger action ap- 

proaches [Sni Eni EH ESI [03] , stochastic quantisation [33,1111122] or lattice gauge 

theory [BSl EZl ESI ESI EOl ED El ESI IB 113 ESI EZl ESI 

troduction see |S2l |83] for textbooks or [HI] for a review about confinement in 
lattice gauge theory). All methods have their merits, and by comparison of 
different results one aims at getting a deeper understanding of the confinement 
mechanism. 

The goal of this thesis is to make another little step further in the understand- 
ing of confinement. A possible mechanism for quark-confinement has recently 
been proposed in [38j, based on DSEs and a power law scaling behaviour of 
Green functions of QCD in the infrared. The authors showed, that this mech- 
anism does not explicitly depend on the Dirac structure of the quark-gluon 
coupling, as Dirac scalar parts yield the same leading scaling behaviour as the 
Dirac vector parts. But the Dirac structure makes the calculation very com- 
plicated. Thus it would be desireable to have a simple model involving bosonic 
particles, that has the same qualitative behaviour in the infrared, i.e. the same 
infrared scaling behaviour of the vertex function. The assumption that this is 
possible is based on the fact, that the string tension is similar for both types 
of particles in the fundamental representation [851 186] . 

In this thesis the system of scalars coupled to the Yang-Mills sector is inves- 
tigated by means of DSEs with respect to its (non-) confining properties, i.e. 
the infrared leading scaling behaviour is determined and compared with the 
scaling exponents in QCD. Herein the potential problem of finite tadpole con- 
tributions at zero temperature will be ignored. The big advantage of this model 
will turn out to be, that despite the fact that the infrared analysis, performed 
here, is complicated by bosonic self-interactions, it should simplify a numeri- 
cal analysis of the mid-momentum region, which could not comprehensively be 
analysed for QCD yet. In addition it would simplify a comparison between the 
Dyson- Schwinger formalism and lattice gauge theory. The point is that on the 
lattice it is very intricate to simulate dynamical fermions, but for bosons this 
can be done more easily. As a result this model may be a suitable test for the 
mechanism of confinement, as it is proposed in [SSI SOI 157] . 

Fundamentally charged scalars coupled to Yang-Mills theory have not been 
dealt with by DSEs yet, but from lattice gauge theory some results for a funda- 
mental Higgs-model in SU(2) are known. From the Osterwalder-Seiler-Fradkin- 
Shenker theorem [SSI ISS] one sees, that there is a free-charge phase, and ad- 
ditionally two other regions, the confinement and the Higgs region, which are 
analytically connected [901 H] • Thus there is no thermodynamic phase tran- 
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sistion [891. I92| 193] . and the distinction must be done by a non-local order 
parameter [OU [95|, |96] , which separates the two regions. 

This thesis will be organised as follows: chapter 2 will give a short introduction 
into the underlying physical concepts, chapter 3 will deal with the specific 
methods that will be used, i.e. DSEs and the infrared power counting analysis. 
The results of my work I will give in chapter 4. In chapter 5 the comparison 
of the scaling behaviour of scalar Yang-Mills theory with quenched QCD will 
be done, before in chapter 6 the conclusion and outlook to future work will be 
given. Conventions, details of the DSEs and parts of the calculation are put in 
various appendices. 
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2. Yang-Mills Theory Including 
Fundamentally Charged Scalars 



In this chapter I will outline the basic concepts of gauge invariance and the con- 
sequences for the Lagrangian density, as the interaction terms and the charge of 
the scalar. I will start with a short introduction into group representations by 
comparing scalar particles in the adjoint representation to scalars in the funda- 
mental representation of the group SU(N). The aspect of renormalisability will 
give further constraints for the possible interaction terms in the Lagrangian 
density. At the end of this chapter I will shortly discuss the model, which will 
be analysed in this thesis. 

2.1. Group Representations of SU(N) 

In the standard model of particle physics quantum chromodynamics (QCD) is 
a gauge theory, i.e. a local quantum field theory which is symmetric under a 
special kind of transformations. For each local gauge symmetry the underlying 
group is a Lie group. For QCD the underlying symmetry group is SU(3). 
Thus the Lagrangian of the system must be symmetric under a so-called gauge 
transformation. 

For one group SU(N) the algebra su{N) can be specified by the commutation 
relation of the abstract generators T" of the group and the structure constants 
jabc q£ ^YiQ group. For the generators of a Lie group the identity 

1 1 ^'^^ 

T^T'' = ^S^hj^ + ^''^')^' (2-1) 

c=l 

holds, where the {f"''"^) d"-^'^ are totally (anti-) symmetric. As the symmetric 
parts drop out in taking the commutator, it leaves the algebra defined by 

[T^ T''] = 2/"''^ (2.2) 

Fields can be described by an (irreducible) representation of SU(N). For SU(N) 
the basic irreducible representation is the fundamental representation. It can 
always be constructed, and further that any other finite-dimensional irreducible 
representation can be generated from it, see [97j for further details. In this 
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representation a field is given by an A^-dimensional vector and the unitary 
matrix representation of a generator T*^ is a x matrix, satisfying the 



commutation relation (2.2). For N > 2 this representation is complex, i.e. 
there is another inequivalent representation for a field ip{x), sometimes called 
antifundamental representation. In order to find group invariants, one can 
always consider products of 

A local gauge transformation can be written in the form 

^{x) -> e'""(")*Xx), (2.3) 

where ip{x) is an arbitrary matter field and a°'{x), a G {1, . . . , A^^ — 1} are 
continuous and differentiable functions, as SU(N) is a compact Lie group. Note 
that a"(a;) is a space-time dependent quantity, i.e. it parametrises a local 
symmetry. 

Another important irreducible representation that exists for all SU(N) is the 
adjoint representation. Given the structure constants f"''"^, which contain the 
entire information about the algebra su{N), one finds that the generators can 
be represented by (A^^ — 1) x [N"^ — 1) matrices given by 

(r")'"= = -iT^^. (2.4) 



2.2. Gauge Invariance 

Forming gauge invariants under an SU(N) of the form iIj{x)iI!{x) is possible 
here, so there are no further problems in constructing mass terms in QCD. But 
a non-trivial theory must have a kinematic part, so necessarily there must be 
derivatives. The problem is, that ordinary derivatives are not gauge invariant 
by this construction, as they act not only on the field but also on the spacetime 
dependent group parameter. Thus they receive additional terms like in 

d^^ix) ^ e'"'^(")*"a^V'(a;) + t (df^a^ix)) t''e'"''(")*Xx). (2.5) 



In order to get a derivative, which transforms as the field itself in eq. (2.3), one 



introduces a so-called covariant derivative, which has an additional space-time 



dependent term, that cancels the redundant term in eq. (2.5), 



d^-^D, = d,-igAlt\ (2.6) 



In eq. (2.6) is the covariant derivative, is the gauge field, and g is the 

coupling strength of the fields ip and Tp to the gauge field. 

Since is covariant, also the commutator [Z)^, Dy\ is covariant. Thus for 
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non-Abelian gauge theories the field strength tensor is defined by, working in 
natural units and Euclidean metric, 

F;, = [D,, DJ\ = d,Al - d,A1 - gr'^A'^Al. (2.7) 

As a result the simplest unitary action, i.e. an action with only two derivatives, 
is 



S= I (2.8) 



Due to gauge invariance we still have to choose a gauge, because the functional 
integral counts all configurations and thus overcounts the ones, that are in the 
same gauge orbit, i.e. which are related by a simple gauge transformation. 
Faddeev and Popov showed a way how to insert this restriction in the funtional 
integral [98] , which I will shortly sketch. 

The gauge fixing is done by a gauge condition (5(/'^(yl)), i.e. considering only 
a hypersurface, such that each orbit subtends the hypersurface only onc^ In 
linear covariant gauges one can express this condition by /""{A) = <9^A^. But 
this insertion changes the measure of the integral. This can be compensated by 
the Faddeev- Popov determinant, which can be rewritten as another contribu- 
tion in the Lagrangiar0 The corresponding particles are called ghosts, which 
are to be seen as scalar particles but underlying fermionic statistics. Note that 
ghosts are given in the adjoint representation, as Faddeev and Popov intro- 
duced them. With this technique of gauge fixing the action of the Yang-Mills 
theory reads 



Sym = I d'x + ^^id^A^Y + ed.Dfc"^ 



(2.9) 



where C, is the gauge parameter, (c") c"" denotes the (anti)-ghost and D'^l' is the 
covariant derivative for ghost fields in the adjoint representatioi]|^ 



D''^ = S^'d^ + gr'^Al. (2.10) 



^In fact problems emerge, if one orbit intersects the hypersurface more than once. These 
"copies" of the same orbit are cahed Gribov copies j99j, thus Gribov's idea was to restrict 
the integration to the first Gribov region, where each orbit only intersects once. Later is 
was shown that still there is overcounting, thus one has to limit the functional integra- 
tion to the fundamental modular region |100) . which causes severe troubles in practical 
calculations. 

^To be brief I only outline this topic, because it can be found in almost every text book on 
quantum field theory, as it is a standard technique. I refer the interested reader to the 
detailed description in e.(;.[101|, as there is also given an analogon to ordinary integration 
techniques, or other didactic presentations fl02l 1103] . 

■^Note that the /"''^ are the generators of the SU(N) of the adjoint representation. 
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2.3. Lagrangian Density for Scalar Yang-Mills 
Theory 

Although pure gauge theory aheady bears enough subtleties, still one has to 
include matter in the Lagrangian to get a theory that is capable of predicting 
processes of the real world. In the standard model of particle physics quarks 
and leptons are the massive elementary particles, whose dynamics and inter- 
actions are to be described. Quarks are known to be fermions, i.e. they obey 
the Fermi-Dirac spin-statistics. In search of a fundamental mechanism of con- 
finement the Dirac-structure may not play the crucial role, but it definitely 
complexifies things, and may conceal the important aspects. Thus it would 
be desireable to have a model, that reproduces the qualitative properties of 
confinement but is free of the complications due to the Dirac structure. 
In QCD the quark-gluon vertex is the link between matter and gauge sector. 
Parametrizing according to Ball and Chiu [IDAl 5 leaves twelve tensor struc- 
tures. Thus it is a system that is very hard to handle. The idea in this thesis 
is using fundamentally charged scalars instead of quarks. Scalar particles are 
bosons, in contrast to quarks, i.e. scalars do not have a Dirac-structure. Look- 
ing at the vertex this means, that there are only two tensor components. So 
the vertex might be easier to handle - first of all in the Dyson-Schwinger for- 
malism, as I will carry out in this thesis, but in lattice gauge theory as well, 
where dynamical fermions are still hard to implement, whereas bosons can be 
simulated efficiently. 

So changing fermions to bosons yields a simpler model compared to QCD in 
the deep infrared region, which is also easier to handle in the mid-momentum 
region. But to be meaningful it furthermore must have the same qualitative 
scaling behaviour. The question is if such a severe modification yields the same 
qualitative behaviouro in the infrared region. 

So the task is to include scalar particles into the Lagrangian. From the physical 
point of view some aspects must be considered: 

• Scalars must be charged, in order to be able to define a conserved charge. 

• Due to the gauge principle the Lagrangian has to be invariant. 

• In what representation does one include the scalars? 

• The Lagrangian must be renormalisable. 

As quarks are in the fundamental representation, also for scalars the funda- 
mental representation is chosen, and from now on I will write the scalar field as 
and its conjugate as 0*. As mentioned above these fields are two inequivalent 
irreducible representations of SU(N). 

So as a result only terms with group invariants are permitted in the Lagrangian. 
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Up to this point there is still an infinite number of possible interaction terms in 
the Lagrangian. But there is another restriction to the possible terms, i.e. the 
theory has to be renormalisable. A quantum field theory naturally generates 
divergencies, as one sums over infinitely many internal modes in loops. But it 
turns out that although divergencies appear, the theory may bring a correct 
physical interpretation about|^ Thus in order to yield a physically meaningful 
theory these divergencies must be cancelled such that one can predict observ- 
able (finite) quantities, at least, after adjusting parameters to the experiment. 
So regularization and renormalisation is a method to render a theory finite 
and predictive. In my opinion the easiest way to understand the idea of renor- 
malisation is by means of so-called counterterms. To get rid of divergencies 
by counterterms, one introduces additional terms into the LagrangiarQ which 
cancel the divergent terms. In order to cancel the divergencies up to all or- 
ders, the counterterms must have the same structure as the divergent functions. 
Another condition for renormalisability is, that there is only a finite amount 
of counterterms. An infinite amount of counterterms would correspond to in- 
finitely many free parameters. Such a theory would not be predictive. Thus 
the renormalisation procedure is not valid any more. 

As a result of renormalisability only several interaction terms in the Lagrangian 
are possible for scalar fields, which can be shown by means of a dimensional 
analysis. Let's start from the fact, that the action 



must be dimensionless (in units of h). The integration gives a length dimension 
of 4, thus the Lagrangian must have length dimension of —4. Due to a Fourier 
transformation one sees that this is equivalent to a momentum dimension of +4. 
With respect to the kinetic terms the dimension of a scalar field must be +1, 
because of the derivative having dimension +1. In contrast to this fermions, as 
e.g'. quarks, have mass dimension |. So a renormalisable theory must only have 
a few distinct interaction terms, each of which have to be gauge invariant and 
renormalisable. These terms are 

• the ones from the Yang-Mills Lagrangian 

• a kinematic one for the fundamentally charged scalars. For gauge invari- 
ance this term must involve two covariant derivatives. These covariant 



^Consider the most intuitive example: the (divergent) self-energy of a point charge in 
classical electrodynamics. But dropping the divergent terms one finds a good description 
of the generated field of the point charge. 

^Note that in fact it is actually a dividing of terms, that are already contained in the 
Lagrangian, rather than adding further terms. 




(2.11) 
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derivatives imply a coupling of scalars to the gauge sector in terms of a 
(two-)scalar-gluon and a two-scalar-two-gluon vertex. 

• a four-scalar coupling with a dimensionless coupling constant, differing 
from the gauge coupling. 

As a first result, the Lagrangian of the scalar Yang-Mills theory reads 

(2.12) 

with the covariant derivatives and the field tensor 

f;, = d,At - d^Al - gr'^A^At. 
For the convenience of the reader I summarise the nomenclature in this work: 

• Greek indices are Lorentz indices, Roman superscripts are color indices. 
Einstein sum convention over identical indices is assumed. 

• F^^ is the field strength tensor. 

• C is the gauge fixing parameter - in this thesis Landau gauge is used, i.e. 

• {(f)*)(f)i is the complex (anti-) scalar field in the fundamental representa- 
tion. Here « G {1, . . . , A^} for SU(N). 

• is the gluon field. 

• (c") c" is the Faddeev-Popov (anti-)ghost. 

• D^^^, Dij^n are the covariant derivatives for the adjoint and fundamental 
representation, respectively, where /"^"^ and t°- are the various generators. 

• m is the mass of the scalar field, and A the coupling of the self-interaction 
of the scalars. 

In the following chapters the generating functional of Green functions will be 
very important for the derivation of Dyson-Schwinger equations. For the full 
Green functions it reads 

Z[J, r, rj, fj, Q] = e-^+^''"^^''^*+''^+'>'+~'^+'^''^^\ (2.13) 
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where (J) J* are the sources of the (anti-) scalars, (77) 77 are the sources of the 
(anti-)ghosts and Q" the sources of the gluons. 
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3. Infrared Behaviour of Scalar 
Yang-Mills Theory 



In this chapter the derivation of Dyson-Schwinger equations is sketched. The 
idea of a skeleton expansion is introduced, which will prove useful in justifying 
the methods, that will be used in later chapters. The relevant DSEs for scalar 
Yang-Mills theory are presented. For later application the power counting 
analysis is defined at the end of this chapter. 



3.1. Derivation of Dyson-Schwinger Equations 

In quantum field theories the quantities of interest are the Green functions, 
also called correlation functions or ra-point functions. Dyson-Schwinger equa- 
tions (DSEs) are the quantum equations of motion for Green functions and as 
they are non-perturbative they serve as a functional method to investigate the 
whole momentum range of QCD. Thus also the infrared region of QCD can be 
analysed, i.e. where the related momenta are far below the momentum scale 
of QCD (Aqcd)- 

F. J. Dyson and J. S. Schwinger separately presented the formalism of DSEs, see 
[31 a ig for the original works or|12ltt3l[Ta[l5in8l[T9liniE2ll3a[Mll371Eg 

for reviews and recent applications in QCD. An alternative method, based 
on equal time commutation relations and Heisenberg's equations of motion, 
is presented in |105j . The underlying idea of Dyson-Schwinger equations is 
that, due to translational invariance, the integral of a total derivative of the 
generating functional vanishes (as boundary terms vanish by definition) 

-^Z[J] = f I?U-^e-^[^l+^'^> = 0. (3.1) 
d(pi J dipi 

Note that Jj is the source of a field yjj, wherein (yj is a reducible representation of 
the irreducible fields and in the index i all color, space-time and representation 
indices are absorbed. S[(p\ is the action of the theory containing ip. From 



eq. (3.1) one can derive a similar expression for the generating functional 



for connected (Vr[J]) or 1-particle-irreducible (F [</)]) Green functions. The 
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modification is based on the relation of tlie generating functional for full Green 
functions to the generating functional for connected Green functions 

Z[J] = e^W (3.2) 

and a Legendre transformation 

r[$] = -W[J] + $,J„ (3.3) 
wherein $j are the averaged fields in the presence of the sources 

= {^^)J■ (3.4) 

The resulting equation 



5S 

6(pi 



dT , , 

+ ^ = (3-5) 



(where A^- is a general propagator of the generic field v?), serves as a starting 
point for the derivation of the DSEs. Note that mixed propagators are con- 
tained until the external sources are set to zero. 

Taking the functional derivative with respect to a scalar, gluon or ghost in 



eq. (3.5) leaves the so-called generating equations for the particles. By taking 
further functional derivatives one obtains the equations of motion for arbitrary 
higher n-point functions. The generating equation must be obtained by doing 
the analytic calculation by hand, but proceeding with this method the deriva- 
tion for higher derivatives becomes a very tedious task. Fortunately there is 
also a graphical method, which I will sketch shortly in the following, for a 
detailed explanation for an algorithm of deriving DSEs I refer the interested 
reader to [106J . where also a Mathematica package DoDSE is presented, which 
derives DSEs by the given algorithm. 

From now on in all diagrams in this thesis a straight line will represent a 
propagating scalar particle, a spring line a gluon and a dashed straight line 
a ghost. Double straight lines stand for the generic field, dots denote bare 
vertices, blobs dressed vertices, crossed circles external fields and squares mixed 



propagators, see table 3.1 
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® 



scalar 



generic field 



gluon 



ghost 



bare vertex 



dressed vertex 



mixed propagator 



external field 



Table 3.1.: Graphical representations for emerging quantities in the derivation of 
DSEs. 



Functional derivatives of the various quantities give 



5$, 
6 



ifc im mjn nki 



r 



]i\...t„ 



(3.6) 
(3.7) 
(3.8) 



These structures can easily by described by graphical means, as can be seen in 



fig. 3J^, where I omitted all indices as the figure serves for pure illustrational 
purposes. The arrows restate the process of taking the functional derivative. 
The last blob stands for an n-point function which receives a further leg. 
After having replaced the structures consecutively one sets the external sources 
equal to zero, i.e. all generic field propagators become gluon, scalar or ghost 
propagators. 
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Figure 3.1.: Replacement rules for the graphical derivation of DSEs. The crossed 
dot represents an external field, which is removed by derivation, a mixed 
propagator, marked as a box, yields a 3-field vertex with two attached 
mixed propagators, and an arbitrary n-point function receives an addi- 
tional leg. 



■ ^ ■ • ■ 

• ^ • 

3.2. Skeleton Expansion 

DSEs are by construction an infinite set of coupled non-linear integral equa- 
tions, as every equation involves also higher n-point functions, which can be 
seen in the diagrams of the following sections. To handle such a system a clas- 
sification is needed, which can be done by means of a skeleton expansion. 
The skeleton expansion is based on an idea of the skeleton of an arbitrary dia- 
gram. Every diagram can be decomposed into its skeleton graph by contract- 
ing all self-energy and primitively divergent vertex contributions to a point. 
Therefore each diagram has a uniquely determined skeleton. But following 
these thoughts one can construct every diagram that can emerge, if one has all 



different skeletons attached with suitable insertions, which are given in fig. 4^ 
for Yang-Mills theory with scalars in the defining representation. 
Applying this method to the DSEs leaves a system of equations, that only con- 
tains different skeletons, each containing one bare and further dressed primi- 
tively divergent vertices. Arbitrary high loop-orders can be obtained by subse- 
quent insertions of the allowed structures. 

The validity of the skeleton expansion is based on some assumptions. First of 
all higher loop-orders must not increase the order of singularity compared to 
the skeleton diagrams. Further the prefactors of the diagrams must "behave 
well", i.e. they do neither vanish nor combine in such a way that different dia- 
grams cancel against each other. Under these assumptions it suffices to analyse 
the first order system via a power counting analysis to gain a self-consistent 
solution about the scaling behaviour of all vertex functions in the deep infrared. 
For illustrative purposes I give a standard example for the skeleton expansion 
of a diagram in the three-gluon vertex, which contains a two-ghost-two-gluon 
interaction, which is not a primitively divergent vertex. Thus it has to be ex- 



panded, and it can be seen in fig. 4.5 that all higher orders can be generated 
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by proper insertions. 

This section shall only give an idea of the skeleton expansion, the explicit ex- 
pansion and further subleties will be defered to the following chapters. 



3.3. Dyson-Schwinger Equations of Fundamentally 
Charged Scalars in the IR-Limit 

As the Yang-Mills system including fundamentally charged scalars shall be used 
as a model for QCD, I want to point out the differences in the basic equations 
already in this section. 

The most obvious difference between quarks and scalars is the interaction of 
these types of particles. Quarks are fermions, whereas scalars obey bosonic 
statistics. Thus the interaction terms in the Lagrangian must be different, as 
explained above, which will influence the DSEs. 

For QCD there are seven primitively divergent vertices, five stemming from 
the pure gauge sector, the quark propagator and the quark-gluon vertex as 
the link of the quark to the Yang-Mills sector. Including scalars instead of 
quarks, two additional vertices come into play, namely the four-scalar and the 
two-scalar-two-gluon vertex, which cause major modifications in the DSEs. All 



nine primitively divergent vertices are listed in table 3.2 The DSE for the 



quark propag atoi0 can be represented graphically, as given in fig. |3.2[ The 



■1 - -1 _ ,@ 



^7^ 
"Si 



Figure 3.2.: Dyson-Schwinger equation for the quark propagator in QCD. 

analogous equation for the scalar propagator is complicated by the additional 
vertices. The scalar propagatoij^DSE in fig. 3^ has eight additional structures, 



each containing a bare two-scalar-two-gluon or four-scalar vertex. Also in the 



gluon DSE additional vertices emerge, as can be seen in fig. 3.4 Nevertheless 



the ghost DSEj^is not affected by the additional interactions up to the graphic 



representation, see fig. 3.5 



The coupling of the scalars to the pure gauge sector is realised by the scalar- 



^See e.g. [7] for the derivation of the quark propagator DSE. 

^For the DSEs in this thesis I derived the various generating equations analyticahy and 
proceeded with the graphical method as explained above. After the publication of |106) 
the DSEs were rechecked with the Mathematica-'packa.ge DoDSE. 

^Note that in QCD the equations for the quark and the ghost are similar, which naturally 
does not hold for the scalar DSE. 
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Table 3.2.: Diagrammatic representation for the nine primitively divergent vertex 
functions of fundamentally charged scalar particles. Coupling scalars 
to the Yang-Mills sector yields 4 additional primitively divergent vertex 
functions. 

gluon propagator vmwmmmwwmfp ghost propagator , 



scalar propagator . ghost-gluon vertex 




scalar-gluon vertex i three-gluon vertex J 



four-gluon vertex ,A four-scalar vertex 




two-scalar-two-gluon 
vertex 




gluon and the two-scalar-two-gluon vertex. The corresponding DSE for the 



scalar-gluon vertex is given in fig. 3.6, the DSE for the latter Green function 
can be found in App. |B} 

Note that only the full Green functions for the propagators are given, all the 
others are truncated in the following way. Due to the complexity of the DSEs 
stemming from the primitively divergent 4-point functions, only the one-loop 
diagrams are considered. The remaining system is skeleton expanded, thus 
there is another diagram in the DSE for the scalar-gluon vertex, involving the 
ghost-triangle and a gluon- "box". So this diagram is of two-loop order. It is 
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Figure 3.4.: Dyson-Schwinger equation for the gluon propagator. 



crucial to take this graph into account, because there is no direct interaction 
between scalars and ghosts in a primitively divergent vertex. Thus there is no 
one-loop diagram that accounts for this interaction type, but in two-loop order, 
such a diagram can occur by interchanging gluons. Thus, leaving this two-loop 
diagram aside would restrict the possible solutions. Also in the two-scalar-two- 
gluon equation this aspect has to be considered. Further details to this will be 
given in chapter |4| 

Already at this point I want to emphasise that this is a suitable truncation. 
The justification for this will be given below. 




Figure 3.5.: Dyson-Schwinger equation for the ghost propagator in QCD as well as 
for a gauge theory including fundamentally charged scalars. 
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3.4. Power Law Behaviour 

Considering the preceding Dyson- Schwinger equations it is obvious that they 
are part of an infinite tower of coupled integral equations. To find a solution 
for the vertices it seems as if one has to solve this coupled system. The authors 
of [I9] showed another possibility to give a statement about the qualitative 
infrared behaviour of the vertices. In their argumentation they pursue the 
following idea. On a classical level Yang-Mills theory is a scale-invariant theory, 
i.e. the action is invariant under a dilatation of the coordinate system and a 
simultaneous transformation of the fields. For a dilatation parameter A the 
coordinate system is rescaled by s — Ax, and the fields transform according to 

3 

Afj_{x) — )■ XAfj_{Xx) for bosonic, and ip^x) — )■ X^tp^Xx) for fermionic fields due 
to dimensional reasons. Therefore the action is invariant 

5* = y (I'^xCym — > J d^xX^CyuiXx) = J dt^x Cym{x') = S. 

The classical scale invariance does not survive the process of quantization, 
i.e. a dimensional parameter is generated, usually called Aqcd, which has the 
dimension of a momentum. The infrared region is the region where the relevant 
momenta are far below the momentum scale Aqcd- As distances are large it 
is a fair assumption due to renormalisation group arguments, that all vertex 
functions should be described as a power law. 

This can be demonstrated by means of a general vertex function F^j^p .., as 



given in eq. (3.9), given by a sum over the different tensor structures T* with 
the proper dressing functions Gi. The p's denote the incoming or outgoing 
momenta in this vertex. 
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For the dressing functions Gi power law ansatze are made as, cf |22j . 

G.ipu ...,Pn)= (4 • • • , 4) mpi ■ ■ ■ ,p^f-'- 
j / 

The dressing functions thus scale as a power 5, of a scaling variable qj , which is a 
function of all external momenta, whereas the index j denotes possibly different 
scaling limits. The prefactors Cij must be constructed in such a way, that all 
vertex functions can be described by the same scaling variable. Note that all 
internal momenta in the integrals must transform into external momenta due 
to dimensional reasons, which can be seen in analytic solutions for two- and 
three-point functions |1U71 11U8] . 

The most obvious scaling behaviour is the uniform limit. Uniform scaling 
describes the fact that for a given Green function there is only one scaling 
variable that vanishes if and only if all external momenta scale to zero uniformly. 



as illustrated in fig. 3/7 This limit is therefore defined as 

o 5 • • • 5 



q^{p\, . . . -> <^ pi, . . . ,pn A tI, • • • , -f constant. 




Figure 3.7.: Uniform limit for an arbitrary vertex function. Plain lines represent 
arbitrary in- or outgoing particles, the blob stands for an arbitrary 
vertex. There is only one vanishing scale in the infrared if all external 
momenta scale to zero. 

In principle the singularities that occur when all external momenta vanish are 
not the only ones that can occur. There can also be divergencies, if only a 
subset of the external momenta scales to zero (so-called soft momenta), while 
others stay finite (referred to as hard momenta). Note that hard momenta are 
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not necessarily smaller than Aqcd to ensure the power law behaviour. 

In the Dyson- Schwinger equations these ansatze make it possible to give a 
statement about the qualitative infrared behaviour of the theory. By insert- 
ing the dressing functions and the tensor structures in the DSEs for the nine 
primitively divergent vertex functions a system of equations for the infrared ex- 
ponents can be extracted. Although an arbitrary DSE contains many different 
tensor structures, it is not the goal of the power counting analysis to find the 
scahng behaviour of the various tensors, but rather to give a statement about 
the scaling of the whole vertex. Thus one aims at finding the infrared exponent 
of the IR dominating tensor, although this analysis does not determine which 
tensor is the leading one. Solving the system of coupled equations for the in- 
frared exponents yields the infrared scaling behaviour of all vertex functions. 
The parametrization of the dressing functions can be done arbitrarily, be- 
cause the scaling exponents are independent of the parametrization. The 
parametrization in this work is done with respect to the bare vertices, which 
can be found in App. |Aj The dressed vertex is parametrised by a product 
of the bare quantity with a dressing function. Thus the system of equations 
involves canonical dimensions, i.e. momentum powers of the bare vertices, and 
anomalous dimensions, that describe the scaling of the dressing function. 
To give some examples for the parametrization, the propagators for the scalar 
particle Sij, gluon D^j^ and ghost D'-^, with the corresponding dressing func- 
tions Bg, Z and G, working in Euclidean momentum space and Landau gauge, 
i.e. taking the limit for the Feynman parameter ^ — )■ in the Lagrangian, are 



shown in eq. (3.10). Similar parametrizations are made for the other primi- 



tively divergent vertex functions. 

- -^."^ . o,M) ^ - ^ , nof,') . -m. 

(3.10) 

The abbrevations for the anomalous dimensions in the uniform limit are chosen 



as given in table 3.3, motivated by the nomenclature in |22|. Further kinematic 
cases, where only a subset of external momenta vanishes are denoted by addi- 
tional superscripts, that symbolise those external momenta that tend to zero, 
whereas the remaining ones stay finite. 
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Table 3.3.: Infrared exponents in the uniform limit of the primitively divergent n- 
point functions. 



n-point function 


anomalous 


n-point function 


anomalous 




dimension 




dimension 


scalar propagator 


Ss 


gluon propagator 


Sg 


ghost propagator 


Sgh 


scalar-gluon vertex 


^sg 


ghost-gluon vertex 


^ggh 


three-gluon vertex 


^3g 


four-gluon vertex 




four-scalar vertex 


§4:3 


two-scalar-two-gluon vertex 


Sssgg 
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In this chapter an infrared power counting analysis for the previously derived 
Dyson- Schwinger equations is done for a massless as well as a massive funda- 
mentally charged scalar particle in Yang-Mills theory. The uniform limit is 
calculated for the full system without any truncation. In the case of a massive 
scalar the importance of kinematic divergencies becomes obvious, whose conse- 
quences for the scaling behaviour of the system are presented in the last section 
of this chapter. This analysis yields self-consistent solutions for the exponents 
of the leading tensor structures in each Dyson-Schwinger equation. 



4.1. Uniform Scaling 

In the Dyson-Schwinger equations for the ra-point functions the powers of the 
particular diagrams are made up by their various constituents, whose powers 
are counted according to the ansatze above. The sum of all contributing expo- 
nents of the ra-point functions and integrals is therefore the infrared exponent 
of the diagram. The easiest way to explain this procedure may be by the help 



of an example, given in figure 4.1 , which is a diagram in the skeleton expansion 



of the scalar-gluon vertex DSE. In figure 4.1 the blobs of the internal dressed 



propagators are drawn explicitly to remind of their contribution in the vertex. 




Figure 4.1.: Ghost contribution to the scalar-gluon vertex. 



The power counting proceeds as follows: Within the diagram in fig. 4.1 there 
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are two loops. This means, that in the integral equation there are two momen- 
tum integrals for the exchanged momentum within the loop. In four dimensions 
one hence gets an additional momentum power of 4 from each loop-integral, 
i.e. an additional exponent of 4. The next contributing parts are the prop- 



agators, and due to their parametrization, see eq. (3.10), they account for 
(^p2^-i+6s ^ ^p2^-i+6g i^p2^-i+Sgh^ rjj_^g vertex powers depend on whether 

the vertex is dressed or bare. By construction of the DSEs there is exactly 
one bare vertex in each term of the equation, and consequently in each dia- 
gram too. The bare vertex has no anomalous dimension but contributes with 
its canonical dimension. With regard to App. |X] the parametrization of the 
exponents is chosen such, that the canonical dimension for the scal ar- gl uon-, 
ghost-gluon- and three-gluon vertex equals | . In the diagram in fig. |4.l| there 
are three ghost-, two gluon- and one scalar- propagator, two scalar-gluon and 
two ghost-gluon vertices. Counting the appropriate parts yields 

((/)') ' {pT'^'' {{pr'^'') ' {{pr'^''") ' (/) ^ ((/) ^^'-^) ' ((/) ^^'^f • 



Regarding only the exponents of the constituents of fig. 4.1 leaves a sum of 
momentum powers 



4 -1+6, + 2{-l+Sg) +3(-l+Sgh) + 



loop scalar gluon 
mtcg. propagator propagators 



;host 

propagators 



bare ghost- 
gluon vertex 



dressed ghost- 
gluon vertices 



scalar-gluon 
vertices 



k + Ss + 25o + 35ah + 25aah + 25. 



In the same way the diagrams in the DSEs in Appendix [B] are counted. 
To be the most singular and thus leading vertex function in the infrared requires 
that the order of the singularity is the highest of all, and thus the exponent of 
interest has to be the smallest of all tensor structures in the DSE. 
Proceeding with the analysis for the theory including fundamental scalar charges, 
already in the power counting analysis three qualitatively different cases have 
to be distinguished the first one is an idealised case that the scalar field has 
no mass initially and no mass is generated either. 

The second possibility is that the scalar particle carries a mass (which in gen- 
eral depends on the momentum). Here another distinction has to be done, 
where after a gauge fixing the residual global symmetry is considered. In the 
Higgs phase the remnant symmetry is spontaneously broken. For a detailed 
investigation of this remnant symmetry see |109j . A suitable quantity for the 

^It is not fair to speak about real "phases" of the system, because according to the Fradkin- 
Shenker-Osterwalder-Seiler theorem [88l [89] there is no thermodynamic phase transition 
between the confinement and Higgs-phases. Thus the word "phases" rather means regions 
in this context. 
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1 I 



2 _ 



Figure 4.2.: Truncated DSE for the scalar propagator. 

distinction of the two phases, wherein the gauge has been fixed to Landau 
gauge [911 [95] is the vacuum expectation value of 



Q= d^x(j)j i^J (fx<P^j . (4.1) 

In the confinement, i.e. unbroken phase {Q) vanishes, whereas it approaches a 
finite value in the Higgs, i.e. broken phase. 

All three phases have different properties, thus each case has to be analysed 
separately. They can be specified by the parameters: 

• massless scalar particle m = (and {Q) = 0) 

• massive scalar particle m 7^ 0, but unbroken phase (Q) = 

• massive scalar particle m 7^ 0, spontaneously broken phase (Q) 7^ 
(Higgs phase) 

In this work the massless and the massive but unbroken cases are studied. For 
these two cases the DSEs are topologically equal, so the diagrammatic repre- 
sentations in the Appendix |B] are valid for both cases. This is no longer true, 
if there is a scalar condensate, because then new primitely divergent vertices 
can emerge. 

The calculation for the systems of equations for the massless and the massive 
scalar are similar for the determination of some infrared exponents. I will 
determine the various exponents for both cases simultaneously and then analyse 
the differences of the two systems in the next subsections separately. 
The diagrammatic form of the DSEs is the same in both cases. Due to a similar 



parametrization, see eq. (3.10 ), of the dressed scalar propagators, only the term 
for the bare propagator is counted differently. Applying the truncation from 
above, the DSE for the scalar propagator simplifies to the terms given in fig. 
1121 

The bare scalar propagator for the massless case is singular, as it scales as 
SW) = ^ (^2)-i_ (4 2) 
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For a finite mass > there is scaling of tlie bare scalar propagator ~ 
as it freezes out at a finite value 

) = 2 I 2 ^ 2 0^ ) ■ (4-3) 

As it will turn out in the calculation it is convenient to define a new quantity 
H, that denotes the canonical dimension of the bare scalar propagator in the 



infrared. The equations (4.2) and (4.3) suggest the definition of fx as 



_ I for a massless scalar: = (a a\ 

^ I 1 for a massive scalar: > 



Using the quantity /i, both cases can be written in one equation (4.5) 



l-6s= min{l - /i, 1 + 6s + Sg + 6sg, 1 + Sg, 1 + 6^} 
-6s = min{-/i, 6s + 6g + 6sg, 6g, 6s} (4.5) 

I want to emphasise that for a massless scalar particle, and as there are no 
other scales than the external momentum, in the coupled system of equations 
for the infrared exponents the canonical dimensions can be subtracted and 
cancel exactly. Thus only anomalous dimensions appear in the solutions for 
the infrared exponents. This is no longer true for the massive case, where the 
canonical dimensions remain in the equations and are involved non-trivially in 
the solutions, which will become obvious in the subsections about the different 
cases. 



Note that due to the parametrization chosen above in eq. (3.10), all the other 
equations (that do not contain the bare scalar propagator) contain the same 
term^ Thus the system of equations for the infrared exponents of scalar Yang- 
Mills theory can be set up by counting the contributing graphs in the truncated 
DSE^ which are given in App. [b} The system is given in App. [Cj 
An important point in one-loop order is to consider the ghost contribution cor- 
rectly. This means, that in the DSEs for the scalar-gluon and the two-scalar- 
two-gluon vertex the diagrams have to be taken into account, that contain a 
two-scalar-two-ghost vertex. As can be seen in the Lagrangian there is no di- 
rect scalar-ghost interaction, but naturally they can interact by interchanging 
gluons. Thus this vertex has to be expanded in primitively divergent vertices. 
Since at one- loop order this vertex does not exist, the leading contributions 
of ghosts have to be diagrams, that interchange gluons between ghosts and 



^In the whole thesis the notations for the infrared exponents for the cases of a massless and 
a massive scalar are equal. Due to the sectional separation it will be clear in the following 
which case is dealed with. 

■^I want to emphasise at this point that the validity of the obtained solutions was also 
checked for the full system of Green functions involving also two-loop diagrams. 
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scalar particles. These two-loop diagrams stem from a skeleton expansion of 
the scalar-gluon scattering kernel. It is important to consider a possible ghost- 
dominance, as it is observed in Yang-Mills theory and quenched QCD. These 
two-loop diagrams must be considered in the truncation scheme. In fig. 4.3 
and fig. the graphs to be included are given explicitly. 





Figure 4.3.: The leading ghost contribution in the DSE for the scalar-gluon vertex. 

This two-loop diagram has to be included into the one-loop truncation, 
to give rise to a possible ghost-dominance. 




Figure 4.4.: Leading graphs of the expansion of the two-scalar-two-ghost interaction 
in the 2-scalar-2-gluon vertex DSE. 

As can be seen in App. [C] the system for the infrared exponents is quite ex- 
tensive. To determine the infrared exponents of the system of equations in 
App. [C| a simplification of the various equations is possible, that is based on 
the assumption of a stable skeleton expansion. The skeleton expansion is an 
instruction how to generate higher loop-order diagrams from first order graphs 



by several replacement rules. Figure 4.5 shows possible graphs in the expan- 



sions of the 3-gluon vertex DSE, that can be produced by insertions given in 



figure 4.6 



According to this the idea that is used to reduce the system is that higher or- 
ders in the skeleton expansion should not become more infrared divergent than 
the first ordei|^ cf. \13\. Otherwise continuous insertions would raise the order 

■*Note that this holds for the assumption, that the prefactors of the diagrams do not vanish 
nor that they combine such that different diagrams cancel exactly. 
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t t 



Figure 4.5.: Skeleton expansion of the 3-gluon vertex DSE. 



of the singularity to an arbitrary high value. Hence the possible insertions into 
first order and consequently higher order graphs yield several constraints, be- 
cause the sum of all contributing exponents must be greater or equal zero. 
Another subtlety is that although higher order diagrams are of the same order 
of divergence as the first order diagrams, the series involving all terms may 
diverge. Thus a convergent series of all contributing diagrams is assumed. 



The insertions for a skeleton expansion are shown in figure 4.6 The constraints 



♦ 



I 




Figure 4.6.: Skeleton expansion elements for given diagrams that raise the loop-order 
to generate arbitrary graphs. 



for the infrared exponents are obtained by counting the IR exponents of the 



extensions in figure |4.6[ The crossed out propagators are already part of the 
initial diagram, i.e. they must not be counted in the constraints. An in- 
sertion raises the order of loops, therefore another exponent representing the 
momentum-integral has to be taken into account in the constraints. 



A power counting of the extension in fig. 4.7 yields 



[P ) [P ) [P ) [P ) ' [P )' [P )^ 
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Figure 4.7.: Exemplary graph from the skeleton expansion extensions. 



Thus the exemplary insertion in fig. |4.7| yields an inequality, which serves as a 
constraint for the infrared exponents: 

-1 + 6s^ -1 + Sg -1 + 6gh + ^ + ^sg + ^ + ^ggh = 



additional ghost propa p-iuon propa- ghost propa- 
loop-integral ^ator 4- ^ ^ 

f o gd.Loi gator gator 



scalar-gluon ghost-gluon 
vertex vertex 



6s + 6a + 6oh + Ssa + Sooh > 



(4.6) 



Counting also the remaining extensions yields nine constraints: 

Sg + 26gh + 26ggh > 

26, + 6g + 26sg > 

25g + Sgh + 5zg + 5ggh > 

35, + 2539 > 

6s + 26g + 6sg + 63g>0 

6g + 6gh + 6sg + 6ggh > 

26s + 6^s>0 

26 g + 54g > 

6s + 6g + 6ssgg > 



(4.7) 
(4.8) 
(4.9) 
(4.10) 
(4.11) 
(4.12) 
(4.13) 
(4.14) 
(4.15) 



There are two further constraints for the value of the exponent for the prop- 
agator of the gluon and the scalar particle, that follow from the bare 4-point 
functions from the theory. 6g > and also 6s > due to a simple consider- 



ation, which is explained by means of fig. 4.8| (note that in all illustrations 



of this kind in this thesis the values that are compared are to be seen as the 
absolute values of the diagrams, because on this level of investigation no state- 
ment about the prefactors is possible.). The infrared exponent of the dressed 
scalar-gluon vertex has to be smaller or equal as the overall infrared exponent 
of the diagram on the r.h.s. in fig. 4.8, which involves a scalar loop, i.e. it 



has a higher or equal order of singularity. Considering the power counting of 
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these two diagrams one obtains the inequahty 6sg < 26s + 6sg- Subtracting the 
dimension of the scalar-gluon vertex leaves the constraint 

Ss > 0, (4.16) 

A similar relation holds for 6g > 0, due to the same arguments, based on the 
fact, that there is a bare four-gluon vertex in the Lagrangian. The infrared 
suppression of the gluon propagator is one necessary condition for the Gribov- 
Zwanziger [110^ 1111] and the Kugo-Ojima |112] confinement scenario. In both 
scenarios the gluon propagator is supposed to be less singular than a simple 
pol^ This is fulfilled by the preceding constraint for all 6g > 0, because the 
full infrared exponent of the gluon is — 1 + 5^. 

Exploiting the preceding constraints in the equations in App. [C| several terms 
drop out, because they are subleading compared to the tree-level term. 

In the tadpole terms, that occur in the equations of the scalar and gluon propa- 
gator the external momentum never enters the inner loop. The tadpole contri- 
butions are thus momentum-independent (infinite) constants that are removed 
in the renormalisation process and cannot be the leading terms compared to 
possible singularities. Therefore, in the actual calculations these terms can be 
dropped in the equations for the infrared exponents of the propagators 6s and 
6g. One obtains: 

-6s = min{-/i, 6s + 6g + 6sg} 

-6g = min{0, 26s + 6sg, 26gh + 6ggh} 

^There is another solution for the system consistent with Sg > 0, the so-called decoupling 
solution, see [HI ES EHl El US], which is observed in lattice calculations [701 [751 [7^ 
and also the refinement of the Gribov-Zwanziger scenario according to [BUlEl]- In this 
case (5g = 1, and all other vertices stay trivial. There is a lively discussion, whether this 
solution is directly comparable to Dyson-Schwinger results, due to complications in the 
definition of the Landau gauge in lattice QCD and in the continuum [ 114| . see [113| for 
a summary of this problem. 



32 



4. Results 



The equation of the ghost is simphfied by renormahsation, cf. [12], [131 [13 [El 
[T7] . which can be done such, that the bare ghost-propagator term drops out. 
Thus the equation for the ghost exponent is uniquely determined as: 

Sgh = -l{5g+6ggH) (4.17) 



The simphfied system is given in eqs (4.18 4.26) 



-6s = min(-/i, 6s + 6g + 6sg) (4-18) 
-6g = mm{0,26s + 6sg,26gh + 6ggh) (4.19) 

-Sgh = li5g + 6gg,) (4.20) 

63g = min(0, 36gh + 26ggh, 36s + '2,6sg, 26s + Sssgg, 25^ + 6sg) (4.21) 

6ggh = mm{0,26g + 6gh + 26ggh) (4.22) 

6sg = min(0, 6s + 26g + 26sg, 26g + 6ssgg, 6s + 26g + 36gh + 26sg + 26ggh) 

(4.23) 

6ig = min(0, 26s + Sssgg, 36g + 6^g + 6ig, 35^ + 6sg + 6ssgg, 36s + 25^g, 

A6g + 3(53,, 45, + 36sg, A6g,, + 35,,^) (4.24) 

= min(0, 26g + 6ssgg, 6s + 25(, + 6sg + ^ss,,, 5s + 26g + 25,,, 

, 25, + 6g + 6sg + 54., 25, + 25, + 35,,) (4.25) 

6ssgg = min(0, 25, + 6ssgg, 36g + 63g + 6ssgg, 5, + 25, + 53,; + 6ssgg, 

, 25, + 6g + 6sg + 6ssgg, 35, + 6sg + 54,, 5, + 25, + 25,,,, 

, 5, + 26g + 6sg + 63g, 5, + 35g + 26sg + 63g, 25, + 25g + 35,g, 

, 25, + 25g + 26sg + 63g, 35, + 6g + 35,,, 

, 25, + 25g + 36gh + 35,g + 26ggh, 



, ^Us I ^^g '-"-'gh ~ '-"-'sg ~ '^'-'ggh: 

, 5, + 36g + 36gh + 26sg + 63g + 26ggh, 
, 5, + 25g + 4:6gh + 26 sg + 36ggh, 26 g + 



)g + A6gh + 36ggh + 5,,gg). 



(4.26) 



This system can be solved by hancj^ when tackled in an effective way, starting 
with the easier parts, and consecutively plugging the determined values into 
the remaining equations. 

For the determination of the infrared exponents I start with the calculation 

^The usual algorithm of computer algebra systems tries all possible combinations. Due to 
the size of this system the computation is very extensive. The system of equations was 
implemented on a desktop computer, but did not yield a result, due to too little memory. 
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of the coupled equations for the ghost-gluon vertex and the ghost propagator. 
From the two non-zero terms in the ghost-gluon vertex equation one drops out 
due to the constraints of the skeleton expansion. Thus 

-^gh = ^i^g + ^ggh) 

Sggh = min( 0, 26g + 6gh + ^dggh) 

Assuming Sggh = 26g + 6gh + 26ggh < to be the IR-leading term shows, that 
this choice is not self-consistent with the ghost propagator equation. If this 
term was the dominating one, the scaling of the ghost-gluon vertex would be 

~^ggh = 25g + 6gh- 

Setting this relation (l.h.s.) equal Sggh = Sg + 26gh from the ghost propagator 
equation (r.h.s) yields 

-2Sg + Sgh = -Sg " 2Sgh =^ Sgh = Sg. 

With this relation each constraint of the skeleton expansion had to be fulfilled. 



Considering (4.7) it is obvious that this constellation is no mathematically 



allowed solution, because 

25g/l + Sg + 26ggh = 26g + 6g + 2{-25 g - 5gh) = SSg < 0, 

that can only be fulfilled hj 6g = 0. This is in contradiction to the initial 
assumption, which was that 26g + 5gh + 25ggh was smaller than zero (and hence 
the dominating part). As a result, the only possible infrared scaling behaviour 
of the ghost-gluon vertex in Landau gauge is the trivial one, i.e. the vertex 
stays bare and scales only with its canonical dimension. This non-enhancement 
of the ghost-gluon vertex, as predicted by [TUl E], can be observed in Yang- 
Mills theory and quenched QCD as well, which was shown both by means of 
DSEs [la [THl EDI EH 122] and lattice calculations |HOl EH [Il5]. Therefore it is 
not relevant up to this point, whether scalars or quarks are coupled to Yang 
-Mills theory, but both particles do not change the scaling behaviour of the 
ghost-gluon vertex, given by 

5ggh = 0. (4.27) 

As an immediate result the ghost exponent is proportional to the gluon expo- 
nent. Applying this relation in the gluon equation leaves 

-5g = min(0, 25s + Sg). (4.28) 
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At this point it is obvious that the gluon equation becomes trivial, i.e. the 
term for Sg is involved in the equation for 6g itself. Thus this terms dominates, 
and no further statement is possible so far, but it leaves a free parameter to 
describe the scaling behaviour of the various Green functions, as it is the case 
in Yang-Mills theory and quenched QCD. In the literature this parameter is 
usually denoted by k and defined by 

^Sg = K= -Sgh, (4.29) 

which is a general relation for linear covariant gauges in Yang-Mills theory 

Hitherto the decisive point of this calculation is that the scalar contributions 
in the equations for the ghost and gluon propagator as well as in the ghost- 
gluon vertex are subleading. Therefore the parameter k, has the same value 
as in Yang-Mills theory or quenched QCD, where the authors of [H] derived 
from the ghost propagator DSE, that this parameter has to be in the interval 
< K < 1, which is in perfect agreement with numerical calculations that 
yield k ~ 0.595 [IHl |35l EI]. Thereby the gluon propagator is suppressed in the 
infrared, whereas the ghost propagator is infrared enhanced. 



The next step is the determination of the infrared exponent of the scalar prop- 
agator 

- 5s = mm{- fi, 5s + Sg + 5sg). (4.30) 

Up to this point there are two possible solutions for 5s, but not both of them 
ensure stability under the skeleton expansion. By a similar consideration as it 
was done for the ghost-gluon vertex the non-trivial value of the scalar propa- 
gator drops out. 

Assume that the second term in eq. (4.30) dominates. Then the comparison 
of the constraint (4.8) 

25s + 5g + 25 sg > 
and the relation for the scalar propagator 

-25s = 5g + 5sg <Q 

would yield an inequality 

-25s > 5g, 

which could only be true iS 5g = 0. But this would be inconsistent with 
the initial assumption —5s<0. As a result the scalar propagator must be 
dominated by its bare value 

5s = fi. (4.31) 
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These results already suffice to determine the infrared exponent of the three- 
gluon and the four-gluon vertex. Inserting the values of the known infrared 
exponents combined with /i < and Sg >0 yields 



Ssg = min( -3k, 3fi + 2Ssg, 2fi + Sssgg) (4.32) 
64g = min( -4k,, Afi + 36sg, 2fi + 6ssgg, 3fi + 6sg + Sssgg)- (4.33) 



Due to the constraint (4.8) there is a lower bound for the terms 



3fx + 2Ssg > fi-2K > -3k (4.34) 
4:fi + 36sg > fi-3K > -Ak (4.35) 



Thus the terms drop out in the appropriate equation. Due to constraint (4.15) 
the similar relation 

2/i + Sssgg > /i - 2k > -3k ( > -Ak ) (4.36) 



holds, so this term cannot be the leading one in both equations (4.32) and 



(4.33). 



At last the constraints (4.8) and (4.15) ensure that 



2/i + 5 ssgg > fl — 2k 



+ =^ 3fl + 5sg + Sssgg > - 3k > -Ak. (4.37) 



I want to point out that this is a very crucial result. At this point all leading 
IR exponents of the vertices from Yang-Mills theory are determined uniquely 
(and independent of the mass of a scalar). Thus coupling scalar particles to 
the Yang-Mills sector does not change the scaling of the pure gauge theory 
vertices, because the scalar contributions do not alter the leading parts in the 
system of equations for the infrared exponents. This is a non-trivial result, 
which could not be expected before this calculation. The persistence of the 
ghost-dominance even with scalar contributions in the vertex functions from 
Yang-Mills theory yields that the result for the leading tensor structure can 
still be obtained by the general result for an arbitrary vertex function, as given 



ineq. (4.38) 



wherein n in the number of ghost-antighost pairs and m is the number of glu- 
ons in the vertex. Analog results have been found in the case of quarks [381 HQ]. 
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Figure 4.9.: (Two-loop) sunset diagram in the scalar propagator equation. 
Already at this point it is possible to determine the general value of 6sg in 



dependence of fx. Inserting all known infrared exponents into eq. (4.23) leaves 



6sg = min(0, 4k + S^sgg, fi + k + 26sg). (4.39) 

With regard to the truncation from above it is not possible to exclude any fur- 
ther terms. But considering that the solutions obtained in this truncation must 
also hold for the full system of equations, a sunset diagram (thus two- loop), 



given in fig. |4.9[ in the scalar propagator equation gives a further constraint. 
Suppose that 4/t + 6ssgg < is the dominating term in the scalar-gluon equa- 
tion. Making a power counting for this diagram, the exponent of the r.h.s. of 
the equation must necessarily be larger or equal to the exponent of the l.h.s. 
to give a consistent solution. Thus 

l + < 4 + 2{-l + 6g) - l + 6s + 6ssgg 

< 4:K + 6ssgg. (4.40) 



The constraint (4.40) is in contradiction to the initial assumption, thus the 
term 4k + Sssgg cannot dominate in the infrared. 

The other two scaling exponents of the scalar-gluon vertex are mathematically 
possible solutions 

Ssg = y -fx-K, (4.41) 
which have also been found in QCD. 



Table |4.1| shows the three mathematically possible solutions of the infrared ex- 
ponents for both, a massless as well as a massive scalar particle. I want to 
point out explicitely at this point, that the graphs containing a four-scalar or a 
two-scalar-two-gluon vertices are subleading in all other Green functions. This 
is an important result of this thesis. The actual values of S^g and Sssgg cannot 
be determined simultaneously for the massless and the massive scalar, because 
there are qualitative differences, which I will explain in the following sections. 
Futhermore I emphasise that this scaling behaviour justifies the truncation pre- 
sented above. By construction all two-loop graphs contain one bare four-point 
vertex. As it has been shown above, the four-scalar and the two-scalar-two- 
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gluon vertices cannot be the dominating terms, and therefore it is fair to neglect 



the two-loop graphs except for the diagrams given in fig. 4.3 and fig. 4.4 Con- 
sidering all two-loop terms in the calculation does not change the final solutions 
of the system, which has been checked explicitly. It is important to note that 
due to this the solutions obtained in the one-loop truncation are not an ap- 
proximation but rather exact self-consistent solutions for the scaling behaviour 
of the leading tensor structures of the full DSEs in the infrared region. 

From this point I will split the calculation of the remaining infrared exponents 
of the four-scalar and two-scalar-two-gluon vertices in two different sections. 
In these sections I will also interprete the physics of the results. 



infrared exponents 










^ggh 


5^g 


hg 


trivial solution 






















partial scaling solution 




2k 


— K 








— 3/« 


-Ak 


full scaling solution 




2k 


— K 


— /i — K 





— 3/« 


-Ak 



Table 4.1.: Uniform infrared fixed points of fundamentally charged scalars. 



4.1.1. Massless Scalar Particle 

In the case of a massless scalar particle the values of and 5ssgg simply follow 
from the lower n-point functions. The four-scalar vertex stays bare up to all 
orders. This can be seen from constraint (4.13). For a vanishing mass of the 
scalar m? = the value of 5^ = (as /i = 0). Thus the constraint yields 



54. > 0, 



(4.42) 



but the value appears in the mm-function of equation (4.25), thus the only 
permitted solution is 

Sis = 0. (4.43) 



The exponent of Sssgg is determined from equation (4.26), which is reduced to 

(4.44) 



Sssgg = mm{2Ssg, Sssgg), 



which yields in combination with the constraint (4.15), that the exponent Sssgg 
depends on the value of Ssg- 

In table 4.2 I complete the solutions for the infrared exponents a massless scalar 
coupled to Yang-Mills theory. 

Mathematically there are three possible solutions for this system, that were 
named according to their scaling properties. I want to emphasise again that 
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infrared exponents 
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full scaling solution 
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— K 


— K 





—3k 


-4k 





-2k 



Table 4.2.: Uniform infrared fixed points of massless fundamentally charged scalars. 



all four-point functions do not dominate in the equations of the lower n-point 
functions, so the scaling behaviour of the leading tensor structures in the Green 
functions is qualitatively the same as for the Green functions in quenched QCD, 
for which the scaling behaviour has been determined in HO] , which gives 
rise to the possibility of a comparison of the two systems on a qualitative level. 
Furthermore it is an interesting result that the scalar propagator is dominated 
by its tree-level term in the infrared. 

The trivial solution stems from the fact that the ghost equation can be fulfilled 
trivially, if the gluon propagator is dominated by its tree-level term. Conse- 
quently K = and this solution cannot be excluded, and hence all exponents 
vanish as the parameter vanishes in the infrared limit. It is not sure whether 
this solution is only a mathematical solution, or if it also has physical content 
for the infrared sector. 

The partial scaling solution describes the Yang-Mills sector as usual according 
to its ghost dominance, wherein the scalars decouple in the IR limit. This 
can be seen in the vanishing exponents of the scalar-gluon and the two-scalar- 
two-gluon vertex, which are the link from the matter to the Yang-Mills sector, 
whereas the gluon propagator is suppressed. 

The full scaling solution yields a scaling scalar-gluon and consistently also a 
scaling two-scalar-two-gluon vertex. The scalar sector couples to the Yang- 
Mills vertices, but this low order of singularity is not high enough to yield a 
linear rising potential. Therefore also in this solution confinement cannot be 
generated. The authors of [38] found that a similar situation is realised in 
QCD. 

4.1.2. Massive Scalar Particle 

The other important case to be presented here is the one of a massive scalar 
particle, but without a scalar condensate, which corresponds to the case of 
unbroken symmetry above. 

This unbroken symmetry ensures, that the Dyson- Schwinger equations as de- 
rived above are not changed by other vertices, that would cause changes in the 
topology of the equations. 

For a scalar particle with a certain mass m?{p'^ — ?► 0) = its propagator 
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freezes out at a finite value in the infrared, in contrast to the divergent bare 
propagator of a massless scalar: 



oc 



2\0 



Thus in this case /x = 1, and therefore (5^ = 1 and 6sg = \/ —1 — k. In 
the following the same nomenclature for all anomalous dimensions as in the 



massless case is used for convenience, cf. Table 3.3 



With the previous results the exponents for the four-scalar and two-scalar-two- 



gluon vertex can be analysed. The equations (4.25) and (4.26) simplify to 



J4s 



ssgg 



6 



ssggi 



min(0, 4/t + 5ssgg, l + 4/t+(0V-l-fi 

, 1 + 4/t + 2(0 V -1 - k), 2 + 2k + (0 V -1 - k) + 
,2 + 4/€ + 3(0V-l-/«)), 

min(0, 3/t + 5ssgg, 2 + 2k + (0 V — 1 — k) + 5ssgg, 

, 3 + (0 V -1 - k) + 54s, 1 + 4/t + 2(0 V -1 - k), 
, 1 + K + (0 V -1 - k), 2 + K + 2(0 V -1 - k), 

, 2 + K + 3(0 V -1 - k), 1 + 2(0 V -1 - k), 5ssgg). 



(4.45) 



(4.46) 



For the rest of the calculation it is convenient to investigate the three different 
solutions separately. 

Trivial Solution 

All other vertices are dominated by their tree-level terms. Inserting the values 



into equations (4.45) and (4.46) yields 



5. 



As 



'ssgg 




0. 



(4.47) 
(4.48) 



as the scalar-gluon vertex does not scale. So also the four-gluon and the two- 
scalar-two-gluon vertex are dominated by their tree-level parts. Also for a 
massive scalar particle it is not clear whether this solution is physically relevant. 



Partial Scaling Solution 

For a scaling Yang-Mills sector, 5s 



1 and 5sg = the equations (4.45) and 



(4.46) become trivial and the exponents are equal to the trivial solution 



-'4s 



'ssgg 




0. 



(4.49) 
(4.50) 
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Therefore this solution describes a decouphng scalar from the scaling Yang- 
Mills sector, as in the massless case. 

Full Scaling Solution 



A naive calculation of the infrared exponents from the equations (4.45) and 



(4.46) is possible, but the solution for the exponent 5ssgg bears an ambiguity, 



as only a range of possible values can be given 



6. 



As 



-1 



K 



'1 -2k < dssgg < 0. 



(4.51) 

(4.52) 



Up to this point three solutions remain, the trivial solution, the partial scaling 
solution and the full scaling solution. The third solution as the only solution 
with a scaling scalar sector suffers from a problem for the four-scalar and the 
two-scalar-two-gluon vertex. The degrees of divergence are small, and consid- 
ering a diagram given in fig. |4.10[ which is part of the fully dressed vertex. 



gives evidence, that the obtained infrared exponent may not be the proper one. 



The simple counting of the uniform exponents of fig. 4.10 shows, that the order 
of singularity is larger compared to the dressed vertex 




> 




Figure 4.10.: Exemplary diagram in dressed four-scalar vertex. 



JAs 



-l + K < -l + 2/t + 2(-l 



Apparently this approximation considering only the uniform limilj^ fails in giv- 
ing the correct values for the infrared exponents 5is and 5ssgg for a massive 
scalar particle. This suggests that soft-singularities may be crucial for the scal- 
ing behaviour, as it will be confirmed below. 



^The calculation was also done with all two-loop diagrams, but did not change the results 
obtained in the one-loop truncation. This is due to the fact that by construction all two- 
loop diagrams contain at least one 4-point interaction, which was shown to be subleading 
in the IR. 
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A similar situation is realised in quenched QCD as the authors of [38] showed 
recently. Diagrams including both soft and hard momenta are the leading di- 
agrams in the IR and kinematic divergencies are related to the mechanism for 
confinement in quenched QCD. The question is whether this mechanism for 
confinement also holds for fundamentally charged scalars in Yang-Mills theory. 
Thus the effect of kinematic divergencies on the scaling behaviour must be 
taken into accounl|£l 



4.2. Kinematic Divergencies 

Hitherto only uniform scaling with its divergences was investigated. Further- 
more there is also another type of singularities that can emerge, referred to as 
soft singularities, as presented in recent publications [221 EH]- The previously 
discussed conformal case is based on the assumption that Green functions can 
be described by a power law in one single scaling variable, which vanishes if 
and only if all external momenta that enter the vertex scale to zero uniformly. 
Thus in the uniform case there is only one relevant momentum scale. 
Soft singularities can emerge, if there are two different momentum scales in- 
volved, but only one of them scales to zero whereas the other one stays finite. 
In the literature vanishing momenta are referred to as soft momenta, whereas 
finite momenta are usually called hard momenta. 

The momentum integral naturally covers all different momentum regions, i.e. 
also finite momenta. Although the momenta stay finite in several kinematic 
regions, it is not guaranteed, that the leading contribution to the loop integral 
arises from soft momenta. In special kinematic cases, where only a subset of 
external momenta vanishes. Green functions may have a different scaling be- 
haviour than in the uniform limit. Note that internal loops can have contribu- 
tions from both scales, therefore considering hard scales requires a modification 
of the former equations for the uniform limit too. 

As an example I will shortly discuss the three-gluon vertex with its soft singu- 
larity as being the clearest one to grasp the idea of soft singularities, due to 
the complete Bose symmetry of the three-gluon vertex. In the uniform case 
all three external gluon momenta vanish, but still the internal momenta can 
be either hard or soft. Due to the hard scales involved, also the equation for 
the leading infrared exponent for the uniform case receives additional terms 
originating from hard-loop contributions. As can be seen in Appendix |X] eq. 



(A. 6), the bare three-gluon vertex reads 



r£.p(p, q, r) = igr'' {6f,,{p - q)p + 5,p{q - r\ + 5p^(r - v)u) <5(p + g + r). 



Note that the authors of [5S] find a different solution based on the formalism of the func- 
tional renormalisation group. 
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Supposing that only one gluon-leg becomes soft yields that the bare vertex is 
dominated by the remaining hard momenta, other possible tensor structures 
that depend only on the soft gluon are subleading in the infrared. Therefore 
for only one vanishing gluon momentum this vertex does not scale canonically 
as in the uniform limit but only with its anomalous dimension. 
A dressed three-gluon vertex with only one soft external momentum possibly 
scales with a different infrared exponent as the uniform exponent discussed 
in the conformal case, thus there is another equation for the exponent ^f^ in 
the soft-gluon limit, that expresses this anomalous dimension for the special 
kinematic case. Due to the symmetry of the three-gluon vertex there is only 
one additional equation, because the scaling behaviour does not depend on 
which gluon is soft. 

The situation becomes more sophisticated for other ra-point functions. In the 
case of the ghost-gluon vertex two additional combinations of soft and hard 
external momenta may reveal different scaling exponents than the uniform ver- 
tex. With either a vanishing gluon- or ghost-momentum there are altogether 
three equations for 6ggh, S^^f^ and S^^f^ to be studied. 

The number of different kinematics rises further for higher n-point functions. 
In general each n-point function can have combinations with n or n — i, i G 
{2, . . . ,n — 1} soft momenta in the infrared. One single hard external momen- 
tum is forbidden due to overall momentum conservation in each diagram, but 
every higher number is allowed. Each kinematic case yields a separate infrared 
scaling exponent for the vertex function. 

For fundamentally charged scalars especially the 4-point functions render a 
huge variety of possible soft singularities, which yield new equations for the 
various IR exponents for each kinematic case. 

Including soft singularities for fundamentally charged scalars yields a system 
of equations, which is far too involved to be handled in this thesis. Fortunately 
this analysis has already been performed in the case of QCD and the similar 
structure of the equations allows to employ these results in the present case. 



As has been shown above scalar contributions in the DSEs of vertex functions of 
the pure gauge sector are subleading i.e. the same structures as for quenched 
QCD, as there are no quark loops, dominate the deep infrared behaviour of 
the pure gauge vertices. The infrared exponents for the kinematic limits of 
Yang-Mills theory are thus unchanged, their values are listed in table 4^, see 
[38] for their derivation. 

With respect to the different kinematic cases the remaining system consists of 



Note that the scalar-loop diagram in the gluon equation is dropped for other reasons. 
Considering this diagram would eliminate any possibility of observing confinement. This 
aspect will be discussed further in the next section. 
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— K, 


2k 





—3k, 


-Ak 








1 - 2k 


1 - 2k 



Table 4.3.: Anomalous dimensions for kinematic cases for pure gauge primitively 
divergent n-point functions. 



the equations for the scalar propagator and the scalar-gluon vertex. In the IR 
limit the contributions to loop graphs in the DSEs from soft loop momenta of 
the order of the small external momentum and hard loop momenta that are 
much larger than these can be separated. The Dyson- Schwinger equation for 



the scalar propagator is given in fig. 4.11 where different possible momentum 
routings in one diagram are given as separate diagrams. The caption s denotes 
a soft momentum, whereas h stands for a finite momentum. 



h 



Figure 4.11.: DSE for the scalar propagator with a diagram for a possible hard 
momentum routed through the loop, s and h stand for soft and hard 
momenta. 



The Dyson- Schwinger equations for the uniform, the soft-gluon and soft-scalar 



limit are given in fig. D.l, D.2 and D.3 in Appendix D 



The determination of the infrared exponents of the vertex in the kinematic 
limits can be done by comparison with the QCD system based on the following 
ideas: 

• The equations of the scalar propagator and the scalar-gluon vertex turn 
out to be topologically equal to the corresponding equations for the quark 
propagator and the quark-gluon vertex in QCD, except for the terms 
including the bare two-scalar-two-gluon vertex. This also holds for the 
four-scalar vertex and the quark- antiquark scattering kernel. 

• The differences in the systems of equations for scalars in Yang-Mills the- 
ory and QCD are based on the canonical dimensions of the vertices. As 



it is explicitly shown in App. [D] the canonical dimensions in the equa- 
tions for the scalars cancel in such a way that the leading structure of the 
equations in QCD is reproduced. In comparison with the calculation in 
the appendix of [3H] the equations for the scalar Yang-Mills theory also 
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yield the same solutions and the values of the exponents S'^ 
anomalous dimension of the scalar-gluon vertex in the uniform limitl^, 
S^g and are the same as of the analogous quantities in QCD. The 
infrared exponents are listed in table [4^ 

• Note that the arguments above hold only because the diagrams containing 
bare two-scalar-two-gluon vertices are subleading, see App. |D] for the 
explicit discussion. 





a" 

sq 


5^ 







-i - K V 


-i - K V 






Table 4.4.: Infrared exponents for the scalar propagator and the scalar-gluon vertex 
in different kinematic limits. 

As a result the statement that soft singularities play a crucial role for the 
infrared scaling behaviour of higher vertices, as it is proposed in QCD, stays 
true for fundamentally charged scalars coupled to Yang-Mills theory. 



1 i 



s 



(a) 
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G * # 

h 



(b) 



Figure 4.12.: Contribution of the four-scalar-gluon vertex to the four-scalar DSE 
with inclusion of kinematic divergences. 



With the strong kinematic divergence in table 4.4 the diagram that dominates 



the four-scalar vertex in the uniform limit is the diagram (b) in fig. 4.12 A 



^^Note the different parametrization than in the calculation in the uniform limit 

for the anomalous dimension of the propagator of a massive scalar. 
^^Due to the difference in the canonical dimensions of the scalar-gluon and the quark-gluon 
vertices the reparametrization q;"^ = S^g + | for the scalar-gluon vertex is necessary to 
map quenched QCD onto scalar Yang-Mills theory. The explicit calculation is given in 
App.lDl 



Note that only diagram (a) of fig. 4.12 has been taken into account in the uniform limit. 
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power counting^ yields 

2 + 2(-l + 2k) + (-^ - k) + 3(-^ - = -2. 

Thus the actual infrared exponent of the four-scalar vertex is much stronger 
than the value obtained in the uniform limit, and gives the stable solution of 



Ms 



-2. (4.53) 



Also the value of Sgsgg is changed by the inclusion of soft-singularities. This 
becomes obvious by the power counting of the dominating diagram, given in 



fie. 4.13 



— * — > r^ m m > • — * 



' (a) " (b) 

Figure 4.13.: Contribution of the two-scalar-three-gluon vertex to the two-scalar- 
two-gluon DSE with inclusion of kinematic divergences. 



2 + 3(-l + 2k) + 2(1 - 3k) + 2(-| - k) = -1 - 2k 

^ dssgg = -1-2k (4.54) 

The physical meaning of these exponents is explained by comparison with QCD, 
which will be discussed in the next chapter. 



-'^•^Notc that the integration over the hard loop does not give a contribution (p^)^, because 
the integral is dominated by momenta that stay finite. 
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In this chapter I will compare the system of fundamentally charged scalars with 
QCD. This is done on a qualitative level, as the infrared scaling behaviour is 
determined by means of a power counting analysis. It will be distinguished be- 
tween the chirally symmetric and the chirally broken phase, as there are severe 
differences in the origin of the scaling behaviour. 



5.1. Chirally Symmetric Phase of Quenched QCD 

It is a suprising and important result of this thesis that the fundamentally 
charged massless scalars feature the same infrared scaling behaviour as chiral 
quarks in quenched QCD. This could not be guaranteed before the calculation, 
as there are fundamental differences in the two systems. As a result of the 
power counting analysis the leading infrared exponents for the primitively di- 
vergent vertex functions in scalar Yang-Mills theory stem from diagrams that 
have analoga in the DSEs in quenched QCD. The solutions that are obtained for 
massless scalar Yang-Mills theory and for quenched QCD in the chirally sym- 



metric phase, are summarised in table 5.1 for the convenience of the reader 



This table shows that the anomalous dimensions of the primitively divergent 
vertex functions are equal, and thus the qualitative infrared scaling behaviour 
of chirally symmetric quenched QCD and the system of massless fundamen- 
tally charged scalars is essentially the same. In this sense quenched QCD can 
be mapped onto fundamentally charged scalars. 



The equality of the systems can be explained by a few arguments. At first sight 
the different nature of interactions seems to be an important difference between 
the theory containing fundamentally charged scalars and QCD. But as it turned 
out in chapter |4] this is only a minor difference in the infrared scaling analysis. 
As already stressed above, scalar particles in the fundamental representation 
are bosons, whereas quarks are naturally fermions. These different kinds of 
particles feature different primitively divergent vertices, and as a consequence 
the system of DSEs for the primitively divergent vertices of fundamentally 
charged scalars coupled to Yang-Mills theory are far more complicated as the 
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corresponding QCD equations, due to the additional four-scalar and two-scalar- 
two-gluon vertices. It was therefore a non-trivial question, if these systems 
had the same infrared power law behaviour. But due to the fact, that these 
additional four-point vertices are subleading, as determined above, the same 
diagrams give the leading contributions, in QCD as well as in the scalar Yang- 
Mills theory. In both systems the pure gauge sector remains unchanged by 
the coupling of further particles. Considering only the anomalous dimensions, 
the exponents of the scalar and quark propagator, respectively the scalar-gluon 
and quark-gluon vertices are equal. 

Remember that due to the four-scalar vertex there is an additional constraint 
5s > for the anomalous dimension of the scalar propagator, see fig. |4.8[ It is a 
crucial fact that there is no primitively divergent four-quark vertex in QCD, so 
that there is no analogous constraint for the quark propagator exponent. This 
has an impact on the solutions of QCD in the chirally symmetric CclSG, clS the 
renormalisation procedure can be done in different ways. Surprisingly, despite 
the restrictive dynamics there are more solutions than for scalar Yang-Mills 
theory. 

trivial solution: In this solution the leading anomalous dimensions of all ver- 
tices vanish. This solution is also found for scalar Yang-Mills theory. 
Although the solution for the scaling exponents seems trivial, it is not 
necessarily a physically irrelevant case. 

scaling solution: Actually this solution contains two different cases, which 
were both found in scalar Yang-Mills theory. In both cases the Yang-Mills 



vertices scale according to the general formula (4.38), but in one case the 
quark sector is decoupled, i.e. the quark propagator and the quark-gluon 
vertex stay bare up to all orders, whereas in the other case the quark- 
gluon vertex is divergent in the infrared, as it scales with {p'^)~'^. So these 
two solutions have the same scaling behaviour as the scalar propagator 
and the scalar-gluon vertex in the partial/full scaling solution for scalar 
Yang-Mills theory. 

other solutions: Due to the complete symmetry of the equations for the quark 
and ghost propagators, the renormalisation can be done in different ways, 
which is the origin for other solutions, where e.g. the quark propagator 
is enhanced instead of the ghost propagator, or both of them. I want to 
emphasise again, that these solutions do not occur in the case of scalar 
Yang-Mills theory. 

Another surprising fact is that there are less solutions in the case of fundamen- 
tally charged massless scalars coupled to the Yang-Mills sector than in chiral 
QCD, although the system "looks" easier for the latter case. I want to point 
out again, that the reason is the additional four-point vertices being subleading 
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in the infrared, and that there is another constraint, which has its origin in the 
four-scalar vertex. 



chirally symmetric 
quenched QCD 








hg 


^ggh 




^qg 




2k 




— K 




—3k 




-Ak 












-K V 


massless scalar 
Yang-Mills theory 


^ 


Sgh 


hg 




^ggh 




Ssg 




2k 




— K 




—3k 




-Ak 












-K vo 



Table 5.1.: Infrared exponents for chirally symmetric quenched QCD and mass- 
less scalar Yang-Mills theory in the uniform limit. 



5.2. Chirally Broken Phase of Quenched QCD 

The feature of confinement can be seen in the potential between two static 
quarks. In the Dyson- Schwinger formalism the relevant quantity is therefore the 
quark-antiquark scattering kerneQ The DSE for the four-quark function can 
be derived according to the algorithm in chapter |3| and it turns out, that this 
equation involves similar vertices as the four-scalar vertex in scalar Yang-Mills 
theor}|^, see fig. B.8 The four-quark DSE is given in [3H], where additionally a 



four-quark-gluon vertex is plotted, which is explicitly given in fig. |5.1| here. As 
this vertex is not a primitively divergent vertex, it has to be expanded by means 
of a skeleton expansion. In the uniform limit the emerging vertices yield the 
same order of divergence for the four-quark vertex as it has been obtained from 
one-loop diagrams. But it turns out that considering kinematic divergencies the 
infrared scaling behaviour differs from the uniform limit, because the scaling 
exponent for the four-quark vertex is raised, as I will sketch shortly. 
In [38] the infrared exponents for the quark propagator and the quark-gluon 
vertex are calculatecQ These anomalous dimensions are summarised in table 



""^Note that there is no primitively divergent four-quark vertex in QCD, whereas there is a 
primitively divergent four-scalar vertex for fundamentally charged scalars. 

^It is necessary to mention that there are more diagrams in the four-scalar DSE in leading 
order, but these graphs involve either one bare four-scalar vertex or one bare two-scalar- 
two-gluon vertex, thus they are subleading in the infrared and they can be neglected in 
the power counting analysis. This justifies the determination of scaling exponents for 
fundamentally charged scalars including kinematic divergencies in Appendix [D] as the 
system can be mapped to quenched QCD. 

•^Here the Dirac scalar parts have the same infrared exponents as the Dirac vector parts of 
the various n-point functions. The mechanism for confinement as it is described in [3 8) 
is therefore not crucially dependent on the Dirac structure of the quarks. This served as 
a motivation for the investigation of fundamentally charged scalars. 
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5.2| in the various kinematic limits, the infrared scahng exponents of the Yang- 



Mills vertices were given in table 4.3 



With these exponents the skeleton expansion of the four-quark-gluon vertex 
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-| -/s: (VO) 





-/« (VO) 



Table 5.2.: Anomalous dimensions for the quark propagator and the quark- gluon 
vertices in the different kinematic limits. Note that the a" for the 
scalar-gluon vertex in the uniform limit was chosen such, that it has the 
same parametrization as S^g here. 

in the limit of hard external momenta, as can be seen in diagram (b) of fig. 



5.i , must be analysed. 




Figure 5.1.: Exemplary diagram in the skeleton expansion of the 5-point function in 
the quark-antiquark Dyson-Schwinger equation, which gives the domi- 
nating infrared exponent. 



The different limits are illustrated in fig. 5^ to clarify the various cases and 
the constituents in the power counting 

In the uniform limit the counting of diagram (a) in fig. 5.2 gives 

3 



(pr Up 



2n4/ /^2\-1+2k 



ip' 



ip' 



(5.1) 



If the momenta of the external quarks are hard, whereas the exchanged gluons 
becomes soft, which means that the quarks are sufficiently spatially separated, 
the power counting gives (for the solution with a scaling quark-gluon vertex) 



(/)^((/)-^^) ((p^) 



4 / X 2 

2\-1+2k 



ip') 



2\~2 



(5.2) 



■*In the uniform limit as described above all external as well as internal momenta are soft. 
^Note that the hard momentum integral does not contribute with a 2 in the exponent. 
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s. s S ^ s 



(a) 



h ^ h ^ h 

i I 

h hSh^ h , h 



h 



(b) 



Figure 5.2.: This diagram emerges from a skeleton expansion of the four-quark- gluon 
vertex in the four-quark DSE. Different hmits yield different scaling ex- 
ponents, "h" denote hard, "s"soft momenta. Diagram (a) shows the 
diagram in the uniform limit, whereas diagram (b) is the leading dia- 
gram in the heavy quark limit with a soft-gluon exchange, which gives 
the correct infrared scaling exponent to yield a confining potential. 



This is a crucial behaviour, because the scaling exponent of the quark-antiquark 
scattering kernel is the relevant quantity for the potential between two heavy 
static quarks. Heuristically confinement means, that only color-singlet states 
can be detected. By separating two quarks of finite mass, the energy between 
the particles will become large enough at some point, such that new mesons 
will be created from the vacuum. In the heavy quark limit this possibility is 
banned. But this implies that the energy between the quarks must be infinite in 
the limit of infinitely separated quarks. Thus a necessary condition for quark 
confinement is a rising static quark potential. For the value of the infrared 
exponent of the quark-antiquark scattering kernel this potential is realised, 
which is shown in eq. (5.3). From the scaling behaviour of the quark-antiquark 
scattering kernel in the heavy quark limit one can find the quark potential V(r) 
by the relation 

/ ipr 
d^P^ ~ |r|, (5.3) 

thus one finds a linearly rising potential for two static quarks. 



As an important result of this thesis, the same mechanism for confinement also 
holds for scalars in the fundamental representation. By means of a power count- 
ing analysis the infrared exponents of the vertex functions in scalar Yang-Mills 
theory can be determined, see App. |Dj where the system of quenched QCD 
is mapped onto scalar Yang-Mills theory. As a result the dominating infrared 
scaling exponents of the vertex functions in scalar Yang-Mills theory are the 
same as for the analogous functions in quenched QCD. Due to this fact one can 
argue that the mechanism for quark confinement also holds for heavy scalars. 
As the DSE for the four-scalar vertex in the heavy scalar limit also involves 
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diagram (b) in fig. 5.2, the dominating scaling exponent for tlie four-scalar 
vertex is also —2. As there is the same scaling exponent there is qualitatively 
the same relation, which was given in eq. (5.3), for the static scalar potential 
as for the heavy quark potential. Through this linearly rising potential heavy 
scalars can be confined by the same mechanism as it is proposed for heavy 
quarks in [38] . 



These results show, that the qualitative scaling behaviour of the infrared domi- 
nating vertex functions of fundamentally charged scalars coupled to Yang-Mills 
theory is the same as for quenched QCD, at least for the cases of massless (mas- 
sive) scalar particles without a scalar condensate, which can be compared to 
the chirally symmetric (broken) phase of quenched QCD, where only in the 
chirally broken phase a confining potential can be seen. Thus this behaviour is 
reproduced for fundamentally charged scalars. 
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6. Conclusion 



In this thesis the DSEs for Yang-Mills theory with a fundamentally charged 
scalar field were derived. Due to the quartic scalar- and the two-scalar-two- 
gluon interaction the diagrammatic form of the system of equations is more 
involved than the system obtained for the primitively divergent Green func- 
tions in QCD. 



For the coupled system of nine primitively divergent vertex functions the power 
counting analysis [19j has been performed. It yields the scaling behaviour of 
the n-point functions in the deep infrared which is qualitatively equivalent to 
that of quenched QCD. 

Considering a massless scalar particle one does not observe confinement. Char- 
acteristic of all mathematical solutions of the system is that (i) the scaling of 



the Yang-Mills vertices in eq. (4.38) is not altered by the coupling of scalar 
particles, (ii) the scalar propagator stays bare and (in) the order of singularity 
of the scalar-gluon vertex is not high enough to trigger a linearly rising po- 
tential between scalars, thus there is no solution that shows confinement. An 
interesting and important reason for this is that the four-point vertices are not 
the dominating parts in the lower ra-point functions. 

For a massive particle there are two different cases, which can be distinguished 
by an either vanishing or finite expectation value of a non-local order param- 



eter Q, given in eq. (4.1). In this work the symmetric case, i.e. {Q) = 
was investigated, where it was shown, that taking into account only the uni- 
form infrared divergences of the vertices is not a sufficient approximation to 
fully describe the infrared behaviour of the system. Assuming that four-point 
functions do not dominate the lower n-point functions as shown in the uni- 
form limit, one finds that considering soft singularities yields a solution with a 
static confinement of massive scalars in the quenched limit. Herein the scalar 



propagator stays bare, and the Yang-Mills vertices scale according to eq. (4.38 ) 



The infrared scaling behaviour of scalars coupled to Yang-Mills theory shows 
that the confinement mechanism as it was proposed by [38l HO], EZ] also holds 
for bosonic particles. This feature suggests that it is not the Dirac structure 
of the quark which is crucial for confinement. It is rather an indication, that 
confinement is exclusively due to non-Abelian gauge dynamics. Thereby, this 
result supports the confinement mechanism in [3Sl SDl 157] . 
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6. Conclusion 



The comparison of the different systems yields that the same diagrams dom- 
inate in the DSEs for the scalar particle, as in the analogous equations in 
quenched QCD. Thus it is reasonable to take fundamentally charged particles 
coupled to Yang-Mills theory as a model for quenched QCD. The advantage of 
this model is that it is easier to find a proper truncation for the DSEs, because 
there are less tensor structures for bosonic particles as for fermions. 
In lattice gauge theories dynamical bosons are easier to implement as dynam- 
ical fermions. Thus an investigation of the presented model brings about a 
severe simplification also in lattice studies, which will be done in the near fu- 
ture |117] . A comparison of the results of both methods will give deeper insight 
in the validity of the mechanism of confinement, as proposed in [38] . 
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■ Bare Vertex-Functions 



The expressions for the bare propagators and vertices can be calculated from the 
quantised Euclidean action of Yang-Mills theory with a coupled fundamentally 
charged scalar field by functional derivations with respect to the in- and out- 
going particles, that contribute to the bare Green functions. This appendix 
gives the bare n-point functions in momentum space, that are obtained by 
a Fourier transformation according to the transformation law for an arbitary 
operator O. 
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A. Bare Vertex-Functions 



bare scalar 
propagator 



So,ij{p^) = -Sij^r^ — ^ 



(A.1) 



bare gluon 
propagator 



(A.2) 



bare ghost 
propagator 



(A.3) 



bare scalar- 
gluon vertex 



rS,ij>(p, q,r) = - -g {t%j r^S{p -q + r) 



(A.4) 



bare ghost- 
gluon vertex 



(A.5) 



bare three- 
gluon vertex 



n%p{P,<l,r) =igr''S{p + q + r) 

- q)p + 5vp{q - r)^ + (5p^(r - p),, (A.6) 



bare four-gluon 
vertex 



^ZipAP,<l,r,s) 



g'^5{p + q + r + s) x 

+rT''%S^,6,, - d^p6,,)\ (A.7) 



bare four-scalar 
vertex 



^o,ijkiiP, q, r, s) = ^5ij5kiS{p -q + r- s) 



(A.8) 



bare two- 
scalar-two- 
gluon vertex 



X Sp.i^S{p - q + r - s) 



(A.9) 
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■ Dyson-Schwinger Equations 



The figures |B]T| [R2l [RSj [RSj (Rej [RTl [RSj |B]9] show only the truncated 
DSEs according to the truncation presented in the main text. Only topologi- 
cally different diagrams have to be taken into account with respect to the power 
counting analysis, because topologically equal ones contribute with the same 
exponents. The "-/-..."in the diagrams denote various permutations, which I 
omit for brevity. Internal propagators are dressed. 

The full DSEs have been checked with DoDSE, a Mathematica-pa.cka.ge for the 
algorithmic derivation of DSES, cf. [106j . 



Figure B.I.: Dyson-Schwinger equation for the scalar propagator. 




Figure B.2.: Dyson-Schwinger equation for the gluon propagator. 




Figure B.3.: Dyson-Schwinger equation for the ghost propagator. 
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B. Dyson- Schwinger Equations 





Figure B.4.: Dyson-Schwinger equation for the 3-gIuon vertex. 



Figure B.5.: Dyson-Schwinger equation for the ghost-gluon vertex. 





Figure B.6.: Dyson-Schwinger equation for the scalar-gluon vertex. 
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Figure B.7.: Dyson-Schwinger equation for the 4-gIuon vertex. 
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B. Dyson- Schwinger Equations 




Figure B.9.: Dyson-Schwinger equation for the 2-scalar-2-gluon vertex. 



60 



m Full System of Equations for the Infrared Exponents of Scalar 
Yang-Mills theory 

The power counting gives a coupled system of equations for the infrared expo- 
nents of the primitively divergent n-point functions. In the uniform limit this 
system is given by: 

1 - 6s = min(l - + 6s + 6g + 6sg, 1 + 6g, 1 + 6s) 
l-6g = mm{l,l + 6g,l + 6s,l + 26g + 63g,l + 26s + 6sg, 

, 1 + 26gh + 6ggh) 
'^-Sgh = min(l, l + 6g + 6gh + 6ggh) 
1/2 + 63g = min(V2, 1/2 + ?,6g + 2^33, 1/2 + ?,6gh + 26ggh, ^2 + 26 g + ^49, 

, 1/2 + 26 g + ^sg, 1/2 + 35, + 26 sg, 1/2 + 26 s + ^ssgg, 

, 1/2 + 2^3 + 6sg) 

^/2 + 6ggh = mmy2,y2 + 6g + 26gh + 26ggh,^/2+26g + 6gh + 26ggh) 

1/2 + 6sg = min(l/2, 1/2 + 6s + 26g + 26sg, 1/2 + 26s + S4s, V2 + 26g + 6ssgg, 
, 1/2 + 26s + Sg + 26sg, 1/2 + 6s + 6g + 6sg, 

, 1/2 + 6s+ 26g + 36gh + 26sg + 26ggh) 

64g = min(0, 26g + 64g, 26s + ^ssgg, 35p + (^Sg + hg, ^6g + 263g, 
, 36s + Ssg + 6ssgg36s + 26sg, 4:6g + S^Sg, A6s + S6sg, 

, A6gh + '36ggh) 

6is = min(0, 26g + 6ssgg, 26s + (^4s, 5s + 26g + 6sg + 6ssgg, 

, 6s + 26g + 26sg, 26s + 5g + ^sg + hs, 26s + 5g + '^^sg, 

, 26s + 26g + 36sg) 

6ssgg = min(0, 26g + 6ssgg, 26s + (^4s, 6s + 6g + 6ssgg, 36g + 63g + 6ssgg, 
, 6s + 26 g + 6sg + 6ssgg, 26s + 5g + ^sg + ^ssgg, ^^s + ^sg + ^is, 

, 6s + 26g + 26sg, 6s + 26g + 6sg + 63g, 26s + Sg + 26sg, 

, 6s + 35,; + 26sg + 63g, 26s + 26g + 36sg, 26s + 26 g + 26 sg + (^sg, 
, 36s + ^9 + 3^53, 26s + 26g + 36gh + 36sg + 26ggh, 

, 6s + 36g + 36gh + 26sg + 63g + 26ggh., 

, 6s + 26g + A6gh + 26sg + 36ggh, 26g + A6gh + 36ggh + 6 ssgg) 
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^) ■ DSEs and Power Counting for the Scalar-Gluon Vertex for 
Different Kinematic Regions 



In this appendix the possibihty of kinematic divergencies is considered. As 
shown in chap. 4.2| several kinematic hmits can exist in a 3-point function, 
which may yield a different scaling behaviour. The exponents for the Yang- 
Mills sector are known and remain unchanged by inclusion of scalars. For the 
scalar-gluon vertex the scaling exponents have to be determined simultaneously 
with the exponent for the scalar propagator. This is done as in the uniform 
case by means of a power counting analysis. 

The Dyson- Schwinger equation for the scalar propagator in the uniform limit, 



which is the only possible limit for a two-point function, is given in fig. 4.11 



The equations for the scalar-gluon vertex in the different kinematic limits are 
given in the figures D.l, D.2 and D.3 Contributions from different momentum 
regions in an integral are expressed in separate diagrams. Captions s and h 
denote soft and hard momenta. 




s ss ss ss ss ss s 




Figure D.l.: DSE for the scalar-gluon vertex in the uniform limit. 



The power counting analysis is performed with respect to the parametrization 
for the scalar propagator and the scalar-gluon vertex as given in the equations 
D.1||D.2 D.3 and D.4 which is chosen with regard to a comparison to quenched 



QCD. Note that the scalar propagator 6'^ in this case is differently parametrised 
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D. Scalar-Gluon Vertex in Different Kinematic Regions 





Figure D.2.: DSE for the scalar-gluon vertex in the soft gluon limit. 




h h , h 




Figure D.3.: DSE for the scalar-gluon vertex in the soft scalar limit, 
as in the calculation in the uniform limit. As it approaches a finite value for 
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D. Scalar-Gluon Vertex in Different Kinematic Regions 



vanishing momenta this parametrization is reasonable. 



S 



sg 



sg 



{p-f 

{pY 



(D.l) 
(D.2) 
(D.3) 
(D.4) 



The usual power counting is combined with the values for the Yang-Mills vertex 



scaling exponents, given in table 4.3 Further a re-definition of the exponent 



of the scalar-gluon vertex in the uniform limit a"^ = 5"^ + 1/2 is done, because 
of the different canonical dimension of the scalar-gluon and quark-gluon ver- 
tices. Thus this definition is in agreement with the parametrization of [38]. 



Altogether this leaves a system of equations, given in eqs. (D.5)-(D.8). This 



system is very similar to the one for quenched QCD, except for the terms that 
contain a bare 4-point function. These terms will be shown to be subleading. 



Si 



a 



sg 
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sg 



(D.5) 



(D.6) 



(D.7) 



mm[0, 2 + 5', + 6g + 5^g, 1 + 61 

mill (1/2; V2 + 4k + S', + 2a^g, 26^^ + 1; 3/2 + 25', + 2k + 2a^g, 
,25tg + l; ,l + 6i + 2K + al, + 1/2; y2 + K + 6i + 2al, 
,2 -2k + 26%) 

mill (O; 1/2 + 4.K + 25%, Z + 5', + 25%, 0; 7/2 + 25', + 25%, 
,1 + 2k + 25%, 0; 1 + 2k + 5%, 5/2 + 6', + 5%, 0; y2 + K + 25%, 
,3 + 5', + 1-2k + 25%, 1-2k) 

5% = mill (0, 3/2 + 2k + 5', + a% + 5%, 3/2 + 2k + 5% + 5%, 5% + 1/2; 
; 1 + 2k + 5; + 5% + a%, 3 + 6', + 25%, 5% + 1/2; 1 + 2k + 5%, 
, 5/2 + 5', + 5%, 0; 3/2 + 5; + 1 + 5% + a%, 5/2 + 5% + 5%, 5% + 1/2 

(D.8) 

As the DSEs in App. [Bjshow there are additional diagrams which do not occur 
in QCD. These diagrams stem from the "sheep" -diagram, which is not part of 
the QCD equations, because of a lack of bare two-quark-two-gluon vertex. The 
graph of interest is given again explicitly in fig. D.4| for illustrational reasons. 
The consequence for the scalar-gluon DSE in the uniform limit is the occurence 
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D. Scalar-Gluon Vertex in Different Kinematic Regions 




Figure D.4.: Additional diagram for tlie scalar-gluon DSE of fundamentally charged 
scalars in Yang-Mills theory, that is not involved in the quark-gluon 
vertex DSE. 



of the following terms 



l + 6i + 2K + al, 



(D.9) 
(D.IO) 



which must be shown to be subleading. The first eq. (D.9) cannot dominate, 



because all terms in the sum of a" are constrained to be > 0. Thus this term 
cannot be leading. 

If the second term were the leading one it would have to be < 0. But for any 



number a < 0: 2a < a. Thus the term in (D.IO) cannot dominate, because 
there is the term 



2^:, + 1 



2 X {5:^ + 1/2) < 5; + 1/2 V < -1/2 



(D.ll) 



which has its analogon in the QCD equation. 

In the soft-gluon limit the terms, that must be shown to be subleading, are 



1 + 2K + 6%, 
5/2 + 5; + 5lg. 



(D.12) 
(D.13) 



The first term (D.12) cannot dominate, because appears in the equation for 
itself in a sum with a positive definite term. Thus it has to be subleading. 



eq. (D.8): 



For the second term (D.13) one can use a constraint from the linear term in 

5/2 + 5', + 5% > 0. 



(D.14) 
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D. Scalar-Gluon Vertex in Different Kinematic Regions 



The two terms under investigation in the soft-scalar hmit are 



l + 2K + d 



(D.15) 
(D.16) 



For (D.15) the idea is equal as in the soft-gluon hmit. S^g occurs in the equation 
for itself in a sum with a positive definite term. So it is necessarily subleading. 



For the term (D.16) there is a constraint from the linear term in eq. (D.7). 
Thus 

1 + 2k + 5%>0. (D.17) 

This analysis implies that the terms in the equations for a"^, and 5^^, that 
have no comparable terms in the QCD equations, must be subleading. As a 
result the system of equations has the same solutions as for quenched QCD, 
namely the infrared exponents of the scalar propagator 6'^ and the scalar-gluon 
vertex in the possible kinematic limits a"^, and 6t„. The solutions are listed 
in the main text in table 14.41 



sg 
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